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Abstract. We give a short, operator-theoretic proof of the asymptotic independence (including a 
first correction term) of the minimal and maximal eigenvalue of the n x n Gaussian Unitary Ensemble 
in the large matrix limit n — » oo. This is done by representing the joint probability distribution of the 
extreme eigenvalues as the Fredholm determinant of an operator matrix that asymptotically becomes 
diagonal. As a corollary we obtain that the correlation of the extreme eigenvalues asymptotically behaves 
like n~ 2 ^ 3 /4p" 2 , where a 2 denotes the variance of the Tracy-Widom distribution. While we conjecture that 
the extreme eigenvalues are asymptotically independent for Wigner random hermitian matrix ensembles 
in general, the actual constant in the asymptotic behavior of the correlation turns out to be specific and 
can thus be used to distinguish the Gaussian Unitary Ensemble statistically from other Wigner ensembles. 

i. Introduction. We consider the n x n Gaussian Unitary Ensemble (GUE) with 
the joint probability distribution of its (unordered) eigenvalues given by 

Vn{h A„) = c n e-%—-% n l A > " A ;l 2 

i<j 

in) (n) 

and denote the induced minimal and maximal eigenvalue by AL? and A mE [ x . Bianchi, 
Debbah and Najim (2008) have recently shown the asymptotic independence of the 
edge-scaled extreme eigenvalues, that is, they proved 

P ( A mm < A mix < v) = P < x ) . p ( Arn t < v) + £»(1) (« - 00) (l) 

with the fluctuations 

A W = 2 l/2 n l/6 ( A W + ^ / = 2 l/2 n l/6 ( A W x _ ^ 

The asymptotic independence can been used (Najim 2009) to design, based on the 
ratio of the extreme eigenvalues, a statistical test for the randomness of matrices that 
does not depend on estimating the actual variance of the distribution of the matrix 
entries (that is, the unknown level of noise in some applications). 

In this paper we shall improve upon these results by showing that the correlation 
of the extreme eigenvalues is a simple, scale-independent device to distinguish the 
GUE statistically from other Wigner random hermitian matrix ensembles (and not 
just from non-random matrices like the ratio-based test). To this end we establish a 
first correction term to the asymptotic independence (1), namely 

P (*& < x. 4t < y) 

= P (*& < *) ■ P (Alt < y) + W-xmy)n- 2 " + 0(n~^) (2) 

as n — > 00, locally uniform in x and y. Here, F2 denotes the Tracy-Widom distribution, 
see (12) below. In fact, the correction term comes as an additional benefit from a short 
and conceptually simple new proof of the asymptotic independence that explains it 
straightforwardly from the asymptotic diagonalization of a certain operator matrix. 
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dimension n 



Fig. i. The dots show the values of the correlation p of the extreme eigenvalues of the n x n GUE as obtained 
from a numerical evaluation of the Fredholm determinant (y) by the method ofBornemann (2008). The dashed line 
shows the leading order term n~ 2//3 /4<7 2 of the asymptotic expansion (3). The circles show the sample correlation 
for 10 6 realizations ofnxn matrices drawn from the GUE. To compare with, the squares show the same for 10 6 
realizations ofnxn hermitian matrices whose algebraic degrees of freedom are uniformly distributed on [— 1, 1], 

In contrast, Bianchi et al. based their original proof on quite a detailed and lengthy 
study of the classical power series of the Fredholm determinants representing the 
probability distributions in (1). 

In Section 2 we discuss the correlation of the extreme eigenvalues. In Section 3 
we comment on universality. Finally in Section 4, we prove the expansion (2). 

2. The Correlation of the Extreme Eigenvalues of GUE. Since both, x) 
and FUy), are probability densities it follows from (2) that the covariance of the 
edge-scaled extreme eigenvalues of the GUE satisfies 

cov (A£L X^ x ) = \n-^ + 0(n~^) (n - co). 

Therefore, because of scale and shift invariance and by recalling (4) below, we get 
the correlation of the unsealed extreme eigenvalues (or of any rescaling thereof) as 

—2/3 

P (AJL A«) = U -^r + 0{n-^) (n -> 00), (3) 

where a 2 = 0.813194792832957- ■ ■ is the variance of the Tracy-Widom distribution. 
Figure 1 visualizes that the leading order term of this expansion is actually quite a 
precise approximation of the correlation even for rather small dimensions n. 

We observe a different asymptotic behavior for random hermitian matrices whose 
algebraic degrees of freedom are uniformly distributed on [—1, 1]. Though the data 
shown in Figure 1 hint at an asymptotic behavior of the correlation of the form 
p ~ cn~ 2//3 here too, the constant c is now, quite distinguishably, about three times 
as large as for the GUE. Therefore, the correlation of the extreme eigenvalues may 
be used as a simple and effective scale-independent device to distinguish the GUE 
statistically from other Wigner ensembles (as defined in Soshnikov 1999). 
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3. Universality. Within the class of Wigner random hermitian matrix ensembles 
there are several limit laws known to hold universally. Examples are the universality 
of the limit eigenvalue density, as given by Wigner 's semicircle law, and of the 
limit distribution of the (properly rescaled) fluctuations of the maximal eigenvalue, 
as given by the Tracy-Widom distribution (see Soshnikov 1999). It is therefore 
reasonable to conjecture the universality of the asymptotic independence of the 
extreme eigenvalues. In fact, the sample correlation (squares) shown in Figure 1 for 
a concrete non-Gaussian Wigner ensemble strongly points into that direction. 

However, since the asymptotic behavior p ~ cn~ 2 ^ 3 observed for this example 
differs in the constant c, it appears that the correction term in (2) has to be specific 
to the GUE. We offer the following explanation for this effect. Choup's (2006, 2008) 
Edgeworth expansion for the GUE, that is, 

P (a£1 < f) = F 2 (f) + 7 (f)«- 2/3 + 0(n" 1 ) (4) 

(where the coefficient y{t) is actually given by an explicit, though quite lengthy ex- 
pression, see Choup 2008, Thm. 1.3), allows us to infer from (2) a likewise Edgeworth 
expansion of the joint probability distribution, namely 

IP (41 < ^max < y) = (1 ~ U-X)) ■ F 2 (y) 

+ ((1 - F 2 (-x)) 7 (y) - y(-x)F 2 (y) + (y)) «~ 2/3 + 0{n^). (5) 

(Note the considerable amount of cancelation that would have taken place within the 
order 0(n~ 2 ^ 5 ) terms if we had established (2) from those Edgeworth expansions at 
the first hand.) Though the leading order term F 2 (f ) of the Edgeworth expansion (4) is 
known to be universal, the coefficient 7(f) of the first correction will in general, as in 
the central limit theorem, depend on some higher order moments of the underlying 
distribution of the matrix entries. Now, since the correction term to the asymptotic 
independence in (2) is contributing to exactly the same level of approximation in the 
expansion (5), namely to the order 0(n~ 2 ^) term, it will also most likely in general 
depend on the specific probability distribution of the matrix entries. 

4. Proof of the Asymptotic Expansion (2). Starting point is the well known 
representation 

1 

p n (A 1/ ...,A„) = -det(JC n (A;,A ; -))^ =1 

of the joint eigenvalue distribution in terms of the finite rank kernel (the second 
equality follows from the Christoffel-Darboux formula) 

k=0 ^ I 

that is built from the L 2 (]R)-orthonormal system of the Hermite functions 
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From this representation we get the determinantal formulae (see Deift 1999, §54) 
P ( X ^ A mil) = det (' " Kn ri2(_oo,X)) ■ P ( A mix < Y) = det (l - K„ ^ (y, 

p f X < Ail, **Y)= det f I - K n r 



min' A max ^5 1 J — uet ^1 - A.„ | l2((-oo,X)U(Y,oo)) 

with the natural constraint X < Y (otherwise the last probability would be zero). 
Whereas Bianchi et al. discuss exactly these determinants using Fredholm's power 
series, we refer to the fact that, for X < Y, 

L 2 ((-oo,X) U (Y,oo)) = L 2 (-oo,X) ® L 2 (Y,oo) 

which implies (for the equivalence of a single integral operator on a union of disjoint 
intervals with a system of integral operators see Gohberg, Goldberg and Krupnik 
(2000), §VI.6.i, or Bornemann (2008), §8.1) 

(») »(») / _ ( 1 ( Kn K " 



P [X < KZ> *£L < Y) = det I I- I £ £ I t"L 2 (-c»,x)ffiL z (Y,oo) I • (7) 



Now, edge-scaling, that is, X = - V2n - 2 _1/2 n~ 1/6 x and Y = \/2n + 2~ 1/2 n~ 1/6 y, 
transforms the kernel entries into 

<>(£, 7) = 2- 1/2 n- 1/6 K„(-V2« - 2" 1/2 n- 1/6 £, -V2n - 2~ 1/2 n- 1/6 )? ) 

JS$ (£,,/) = 2- 1/2 n- 1/6 K„(-72n - 2- 1/2 «- 1/6 £, v^n + 2- 1/2 n" 1/6 )7 ) 

4? (£17) = 2~ 1/2 n- 1/6 K n (V2n + 2- 1/2 n- 1/6 Z, -Vln - 2' 1/2 n' 1/6 r,) 

K&\s,tj) = 2- 1/2 n- 1/6 K„(V2n +2- 1/2 n- 1/6 ^, V2n + 2- 1/2 n- 1/6 ^) 

and we obtain (note that X < Y eventually for n large, if x and y stay bounded) 

P (" A min < *) = det (' " ^ K 1? P *) / P ( A mix < y) = det (l - P y K$ P y ) , 

p^o.a^)^,^'-* ( 8 , 

Here, the operator matrix operates on the space L 2 (R) ® L 2 (R) and the orthonormal 
projection Pf : L 2 (]R) — ► L 2 {t, 00) is simply given by the multiplication operator with 
the characteristic function Xt of (f, 00). Plugging the Plancherel-Rotach expansion 
(Szego 1975, Theorem 8.22.9) of the Hermite functions, that is, the locally uniform 
expansion 

<p n (Vln + 2- 1/2 n- 1/6 f) = 2 1/4 «- 1/12 (Ai(f) - \Ai'(t) n - 1/3 + 0(«- 2/3 )) , 

into the rightmost expression defining the kernel K n in (6) yields the locally uniform 
asymptotic expansion 

K^&r,) = K&\s,tj) = K& V ) + 0(n" 2 / 3 ) (n - 00) 
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with the Airy kernel 

Furthermore, Choup (2006, Theorem 1.2) proved that this expansion of kernels 
implies the asymptotic expansion 

P t K^Pt = P t K$P t = P t KP t + 0(n- 2/3 ) (n -> 00) 

of the induced trace class operators (that is, the error 0(n~ 2 ^ 5 ) is also valid in 
trace norm). Completely analogously, by using the Plancherel-Rotach expansion 
once more and recalling the symmetry <p„(—t) = ( — l) n (p n (t), we obtain the locally 
uniform expansion 

K[f(C,r,) = K^\c,v) = H- 1 ) n ^ Ai ^) Ai (v)n- 1/3 + 0(n- 1 ) (n - 00), 

which, by the same arguments as Choup's, implies the validity of the asymptotic 
expansion 

P t K[l ) P s = P t Kffp s = l(-l) n - 1 (xtAi®Ai X s)n~ 1/3 + 0(n- 1 ) (n 00) 

of the induced trace class operators. This shows that the off-diagonal operators in (8) 
have trace norm 0(n~ 1/ ' 3 ). Thus, by the local Lipschitz continuity of the determinant 
with respect to the trace norm (see Simon 2005, Thm. 3.4) we get 



P 



= det (l - P, K$P X ) ■ det (i - P y K^P y ) + 0(n" 1/3 ) 

= P(-S&<*).p(A&<y)+0(«- 1 /^ 

where we have used the multiplication rule of the determinant for diagonal operator 
matrices. This proves the asymptotic independence result (1) of Bianchi et al.; here 
with an additional estimate of the order of approximation, however. 

We now go one step further and make the error term explicit to the leading order. 
We start with the factorization 



P 



(-A^*,A max , y )=det^ q 



det 2 



(I-P x K^P x )^P x K^P y s 
\-1d 



As above, the first determinant evaluates to 

det (i - P, Kg>P x ) ■ det (2 - P y X^P y ) = P < *) • P (A^ix < y) . (9) 
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The second determinant can be written for short as 

dc Jj (-1)^-1/3 / T xy + 0(n-^)\\ 

{ 2 { Tyx + 0(n-W) ) ) 

with the rank-one operator 

T u = (I - PtKPt^XtAi <g> Ai Xs - 

By a straightforward operator decomposition (see Gohberg et al. 2000, (1. 3. 7)), this 
evaluates and expands to 



det (i - \n- 2/3 (T xv + 0(n- 2/3 )) (t ijx + 0(n 



-2/3- 



det (l - ln- 2/3 T xy T yx ) + 0{n 



4/3^ 



The last determinant can actually be evaluated exactly. Indeed, by using the fact that 

(/ ® S) ( v ® w ) — (g> v )f ® w > m d thus 

det(I- (f®g)(v®w)) =det(I- (g,v)f®w) = 1- (g,v)(f,w), 
we obtain 

det (j - in- 2/3 T,. y T yi ) = 1 - \u{x)u{y) n~ 2/3 

with the function 

u(t) = {(I - PtKPt^xtAUAixt)- (10) 
To summarize, our result so far is 

p (-41 < *. ASL y) 



(41 < y) (1 - 1« W«(y) ^ 2/3 + o(«^ 3 )) . 



The product of the probabilities with the term u(x)u(y) can further be simplified by 
expanding the first factor, written as the determinantal expression in (9), through 

det (i - PtK^Pt) = det (i - PtK&Pt) = det (I - P t KP t ) + 0(n" 2/3 ). (11) 



Now, by introducing the Tracy-Widom (1994) distribution 

F 2 (f) = det(J - P f KP f ), (12) 

and recalling that the function u(t) defined in (10) actually satisfies w(f) = F 2 (f) / ^(f) 
(see also Widom 2004, p. 1132) — a formula that was obtained in course of Tracy and 
Widom's derivation of their famous Painleve II representation of F2 — we finally get 

p < */ A mix < y 

= p < *) ■ p < y) - i«(y) «" 2/3 + o(n- 4 / 3 ) 

which is easily seen to be equivalent to the asserted expansion (2). 
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Remark. Numerical experiments using the methods of Bornemann (2008) show 
that the error term of this expansion is indeed not better than of the order 0(n~ 4/3 ). 

Acknowledgements. The author thanks Herbert Spohn for helpful comments 
on a first draft and for bringing up the issue of universality. 
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